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Abstract 

We construct the gravitational dual, in the Unruh state, of the "jammed" phase of a CFT at 



dn- 
strong coupling and infinite N on a fixed five-dimensional rotating Myers-Perry black hole with 

equal angular momenta. When the angular momenta are all zero, the solution corresponds to the 

five-dimensional generalization of the solution first studied in [1]. In the extremal limit, when the 

angular momenta of the Myers-Perry black hole are maximum, the Unruh, Boulware and Hartle- 

Hawking states degenerate. We give a detailed analysis of the corresponding holographic stress 



energy tensor for all values of the angular momenta, finding it to be regular at the horizon in all 
cases. We compare our results with existent literature on thermal states of free field theories on 
black hole backgrounds. 
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I. INTRODUCTION 

Since the revolutionary discovery that black holes radiate with an almost thermal spec- 
trum, much effort has been devoted to gaining a deeper understanding of this process. While 
considerable advances have been made since Hawking's seminal paper 2] , there are still am- 
ple open issues that remain to be addressed. In particular, most of the work regarding 
Hawking radiation has focused on free fields; studying the behavior of interacting fields 



i some progress has been made with weakly 
3|, strongly interacting fields pose a much 



around a black hole is more challenging. Thoug 
interacting fields using perturbative methods 
more difficult problem. 

Fortunately, the AdS/CFT correspondence 4| provides an invaluable tool for tackling this 



problem. Early applications of holography to studying Hawking radiation coupled the field 
theory to dynamical gravity, which is holographically dual to Randall- Sundrum braneworld 
models [5H7| ■ However, if one is not interested in the backreaction of the Hawking radiation 
on the spacetime, the CFT can be placed on a nondynamical background spacetime DM.. 
Then AdS/CFT claims that the CFT is dual to a spacetime Ai with conformal bound- 
ary dM., with M. a solution to classical Einstein gravity with negative cosmological constant. 
While most applications of the duality place the field theory on ordinary Minkowski space, 
one may in principle select any conformal boundary structure one desires. In particular, in 
order to study Hawking radiation in the CFT, [8|, |9| considered a CFT living on fixed d- 
dimensional background containing a black hole of size R and temperature Tbh- Far away 
from the black hole, the CFT is in a thermal state with a prescribed temperature T^. From 
the perspective of the CFT, the black hole and the plasma at infinity act as heat sources 
and sinks, and one can study the exchange of heat between them. 

Regarding the dual geometry, [8|, |9[ conjecture two possible families of solutions: so-called 
"black droplets" and "black funnels," illustrated in Figure [TJ In the black funnel solutions, 
the horizon of the boundary black hole is connected to an asymptotically planar horizon 
in the bulk; this connectedness manifests itself in the CFT as strong coupling between the 
boundary black hole and the heat bath at infinity, leading to an exchange of heat between 
them. In the black droplet solutions, on the other hand, the horizon of the boundary 
black hole is disconnected from the asymptotically planar horizon in the bulk, implying 
that the boundary black hole is not coupled to the plasma at infinity (or rather, that the 
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FIG. 1: A sketch of the relevant solutions: (a): black funnel and (b): black droplet above a 
deformed planar black hole. Both describe possible states of dual field theories in contact with 
heat baths at temperature Too on spacetimes containing black holes of horizon size R. The top line 
corresponds to the boundary, with the dots denoting the horizon of the boundary black hole. The 
shaded regions are those inside the bulk horizons. 



coupling is suppressed by 0(1/N 2 ) is this large-iV picture). From the perspective of the 
CFT, the transition from funnels to droplets is reminiscent of many soft condensed matter 
systems which exhibit a transition from a fluid-like behavior to rigid behavior with no 
flow (e.g. sand in an hourglass, cars on a highway). To borrow from the soft condensed 
matter nomenclature, we will refer to the the CFT transition from funnels to droplets as a 
"jamming" phase transition, and will denote the CFT dual of the droplet as a "jammed" 
phase 1 . 

The authors of fl 3 postdate that this januning transition mig ht occur a 8 the s 1Z e * 
of the boundary black hole (or alternatively, as the temperature T^ of the heat bath at 
infinity) is varied. Thus the dimensionless parameter RToo should characterize which phase 
is thermodynamically preferred. 

An analytic construction of black droplets and black funnels has only been performed in 
select few cases. In 10|, 111] , black droplets and funnels were constructed from the AdS C- 



metric, while 12J constructed analytic funnels in d = 2 dual to the Unruh state of the CFT 
(black droplets cannot exist in d — 2). However, many of the most interesting cases (notably, 
asymptotically flat boundary black holes in d > 2) do not lend themselves to analytic 



1 We thank Jean Carlson for introducing us to this terminology. 



construction, and one must resort to numerics. To that end, [l3|] numerically constructed 
black funnels dual to the Hartle-Hawking state of the CFT, while |l4| constructed so-called 
global funnels to study the exchange of heat between two boundary black holes of different 
temperatures. 

Slightly less effort has been made to construct black droplets, however. Besides those 
obtained from the AdS C-metric, the only construction to date has been that of [1], which 
constructed a droplet dual to the Unruh state of a CFT on a d = 4 Schwarzschild background 
(hereafter generally denoted Schw^). In our notation, this solution sets T^ = 0, so that the 
bulk geometry is that of a black droplet suspended above an extremal Poincare horizon deep 
in the bulk. The boundary stress tensor contains no flux terms, and the nonzero components 
of the stress tensor exhibit a 1/r 5 power-law falloff at large distances from the black hole. 
This relatively rapid falloff is characteristic of the jammed phase, as it indicates that the 
CFT plasma is relatively well-localized near the black hole. Interestingly, the stress tensor 
is regular on both the past and future horizons; the authors speculate that the inclusion of 
one-loop graviton corrections in the bulk would render the stress tensor singular on the past 
horizon, as is typical of Unruh states. 

Our purpose in this paper is to generalize the result of [if by adding a new ingredient: 

nn 

rotation. Though [9I.I15I| compute the stress tensor on a rotating BTZ black hole background, 
the effect of rotation on a stress tensor in an asympototically flat black hole spacetime has 
not been studied extensively. One might naturally expect that the inclusion of rotation 
can be accomplished by generalizing the result of Jj to a Kerr spacetime, but because the 
Kerr metric is cohomogeneity two, this would lead to a cohomogeneity three problem in the 
bulk. Instead, we take the boundary spacetime to be the d = 5 equal-angular-momentum 
Myers-Perry metric 16[ (see also [17| for a review), which is known to be cohomogeneity 
one and thus leads to a cohomogeneity two problem in the bulk. 

The generalization to a non-extremal rotating black droplet should not be expected to 
exhibit very different behavior from its non- rotating counterpart: the boundary black hole 
still acts as a heat source at finite and nonzero temperature, and one studies the thermal 
coupling between the black hole and a zero-temperature heat bath at infinity. However, the 
presence of rotation provides us with a new parameter to tune, which we can use to take 
the black hole to extremality. The fact that both the extremal black hole and the heat bath 
at infinity are at zero temperature implies that, in some sense, the CFT is in its ground 



state, and thus can be thought of as being in the Boulware vacuum. Similarly, since in the 
extremal case Tbh = T^, one can think of the black hole and heat bath at infinity as being 
in thermal equilibrium, and thus the CFT can be thought of as being in a Hartle-Hawking 
state. Thus the extremal droplet is dual to a degenerate state of the CFT corresponding to 
all of the Unruh, Hartle-Hawking, and Boulware states. 

This paper is organized as follows. In Section [Til we review the behavior of stress tensors 
on black hole spacetimes, and calculate the form our stress tensor must take on a Schws 
background (i.e. when the angular momentum of the Myers-Perry black hole is taken to 
zero). In Section HTT| we briefly review the numerical method used and discuss the numerical 
construction of our solution. In Section HVl we compare our numerically extracted stress 
tensor to the expectations in Section [TTJ and discuss future directions to pursue. 



Note: In the final stages of this work we learned of [18j|, which also constructs rotat- 
ing black droplets and may have some overlap with our work. Their paper will appear 
simultaneously with ours on the arXiv. 

II. A REVIEW OF STRESS TENSORS IN CURVED SPACETIME 

A common difficulty of field theories in curved spacetime is the lack of an unambiguous 
notion of "particle." This poses no problem, however, if one instead limits oneself only to 



currents. As discussed in [191 ]. when a physical system is probed, the objects that couple to 
probe fields are currents. In other words, the currents are the ojbects that appear in the 
interaction terms in the equations of motion of an interacting field. In addition, a current 
is a genuinely local object and should retain a well-defined meaning even in the presence 
of strong external fields (such as spacetime curvature). Thus most of the work devoted to 
understanding the process of radiation from black holes has focused not on quantum fields 
themselves, but on the renormalized vacuum expectation value (v.e.v.) of their stress-energy 
tensor (T^"). 

In particular, the stress tensor of quantum fields can exhibit properties that are forbidden 
of classical stress tensors. For instance, though a classical energy density must be non- 
negative, it is well-known that even in flat spacetime, the local energy density of a quantum 



field need not obey this same restriction 2 . Perhaps the most famous example of a negative 



energy density is the Casimir effect 



20] , though the existence of negative energy densities 



in local quantum field theory is in fact quite general 21]. 



In the context of general relativity, 19] studied the behavior of the energy-momentum 
tensor of a massless scalar field in a two-dimensional spacetime containing an accelerating 
mirror. While the total energy radiated by the mirror was positive, a negative energy density 
and energy flux were present in the region near the mirror, whereas far from the mirror 
the stress tensor took the classical form of an outward flux of positive energy. A similar 
result was found by [22], which calculated the renormalized v.e.v. of the stress tensor of a 
massless scalar field in a two-dimensional model of black hole collapse. As in the case of the 
radiating mirror, the near-horizon region of the newly-formed black hole was characterized 
by a negative energy density and a flux of negative energy flowing into the black hole, while 
far from the black hole the stress tensor again took the form of an outward flux of positive 
energy. These results are consistent with Hawking's picture of particle-antiparticle creation 
near the horizon of a black hole, where negative-energy particles flow into the black hole 
while their positive-energy partners radiate to infinity: by conservation of energy, a positive- 
energy flux at infinity should be accompanied by a negative-energy flux through the black 
hole horizon. 

These exact results were obtained by exploiting the conformal flatness of two-dimensional 
spacetimes. Nevertheless, good approximations to thermal stress tensors have been obtained 
in higher dimensions. Notably, [23] used the Bekenstein-Parker Gaussian path integral ap- 
proximation [24| to find the approximate form of the stress tensor of a conformally invariant 
scalar field in the Hartle-Hawking state on the (four- dimensional) Schwarzschild background. 
The behavior of the stress tensor was similar to the two-dimensional case: at infinity, the 
stress tensor behaved like a classical thermal stress tensor, while sufficiently near the hori- 
zon (for r < 2.3AM) the energy density became negative. These results were consistent 
with earlier numerical calculations [25[ of the stress tensor on the bifurcation two-sphere of 
Schwarzschild. However, while one naturally expects a negative energy flux near the horizon 



2 Here, we define the local energy density as measured by an observer with velocity u^ as 

p = vPu v {T, a ,) 



in an Unruh state, it is not clear whether the ubiquity of negative energy densities near the 
horizon should extend to Hartle-Hawking states as well. 

All of these results focused explicitly on conformally invariant noninteracting massless 
scalar fields. One might therefore imagine that they have little bearing on the stress tensor 
we will obtain for our strongly coupled CFT. In fact, the jammed phase dual to the droplet 
we construct cannot exist in a free field theory, so a priori none of the results listed above can 
even approximately predict the behavior of the stress tensor we will extract. Nevertheless, 
it is possible (and indeed, quite probable) that the stress tensor of quantum fields on black 
hole backgrounds should exhibit some kind of universal qualitative behavior. This universal 
structure was studied in detail by Christensen and Fulling [26|, who sought the general 
behavior of a static, spherically symmetric stress tensor on Schw4 by solving the conservation 
equations 

V M T^ = (2.1) 

directly, without restriction to any particular kind of field theory. Their results are partic- 
ularly useful to our case, because they are also valid in the case of a strongly coupled CFT. 
In fact, Schws is just the zero-angular momentum limit of the Myers- Perry black holes that 
we Schwarzschild in higher dimensions. 

A. The General Static Spherically Symmetric Stress Tensor on Schw^ 

Consider the d- dimensional Schwarzschild metric Schw^ (also called the Tangherlini met- 



dsl hw = -f(r)dt* + ^- + r*dn 2 d _ 2 , /(r)^l-(^)"" 3 . (2.2) 



ric, after its discoverer |27l|): 

f(r) 

This spacetime is qualitatively identical to the familiar four-dimensional Schwarzschild, with 
horizon located at r = r$ and a temperature T = /'(ro)/47r — (d — 3)/47rro. Now, a static, 
spherically symmetric stress tensor must take the form 



T% 



(t\ t\ ^ 

T r t T\ 

v o o rV;/ 



(2.3) 



where all components are functions of r only and the indices i, j run over the angular coor- 
dinates. Inserting this form into the conservation equations (j2.ip and using the metric (I2.2I) . 



we obtain the following differential equations: 

o = a r r, + ^r„ (2.4a) 

= d ,.T* r + (*Zl - d -f^ r °lfr 2 ) T% - d ^Kl^tl n _ izl r%. (2.4b) 



r 2r f(r) J 2r /(r) ' r 

Equation ( j2.4a|) can be immediately integrated to yield 

T r t = K^Y~ 2 . (2.5) 

To integrate equation fl2.4bl) . we first substitute T l t = T M M — T r r — (d — 2)T%, and obtain 



f{r) 



1 /' r <7r' 

- / ((d - 3)T%(r') + (d - 2) (2(r'/r )^ 3 - d + f ) 1^(0) — 



• (2-6) 



Equations ( I2.5P and (I2.6P express the most general static, spherically symmetric stress tensor 
on Schwrf in terms of two arbitrary constants Q and K and two arbitrary functions T M M 
and T u q. In order to study the physical behavior of these solutions, it will prove useful to 
define 

0(r) = T%(r) - — 1-_ ^(r), (2.7a) 

d — 2 r dr' 

G(r) = -—- / {2(r'/r ) d - 3 - d + l) 9(r') — , (2.7b) 

^ Jr r 

H(r) = \f (^ + (r'/r o y-*) T%{r>) d —. (2.7c) 



2j ro \ 2 ^ ' u/ / ^ ' r 

Converting to the tortoise coordinate dr* = dr / fir), the stress tensor can then be expressed 
as the sum of four pieces: 

t\ = (T (1) r u + (T {2) r u + {T m y v + (T {4) r v , (2.8) 



where (in t, r* coordinates) 

W, = <N! {-^P if (r) + i n(r), ^P ff(r), ^ 7v» *<} 



PW\ = k (r " M 



d-2 



f(r) 



(l 1 o\ 

-1 -1 
\ 0/ 



(2.9a) 
(2.9b) 



(T( 3 ))^ = diag {-^P G(r) - (d - 2)0(r), ^^ G(r), 0(r) 5}| , (2.9c) 

( T (4)) M , = (r ° / / [ r ) ) ' 2 diag {-1, 1, 0} . (2.9d) 

By converting to ingoing Eddington-Finkelstein coordinates dv = dt + dr*, it is straight- 
forward to show that the stress tensor is regular on the future horizon only if Q vanishes 3 . 
Thus the above decomposition is convenient because each of the T/{\ isolates some interesting 
physical behavior of the stress tensor: Tm is the only of the Tu\ with a nonzero trace (and 
in fact, it is only a function of the trace T^); Ti 2 ) contains nonzero flux terms T r * t ] ^(3) is 
the only of the Tu) which is both traceless and has nonzero tangential pressure components; 
and T(4) is singular on the future horizon. 

Christensen and Fulling proceed to make use of these results to study the behavior of the 
stress tensor of fields in the Unruh, Hartle-Hawking, and Boulware vacua on Schw4. Our 
goal, however, is less broad: we only wish to gain some insight on the v.e.v. of the stress 
tensor of the jammed phase of the CFT. To that end, we make the following observations 
regarding our droplet: 

• The trace anomaly of any odd- dimensional CFT vanishes. Since we work in d — 5, we 
expect T(i) to make no contribution to our stress tensor. 

• According to the arguments reviewed briefly at the end of Section [J we expect black 
droplets to be dual to a jammed phase of the CFT. Such phases are characterized 
by a suppressed exchange of heat between the black hole and the thermal plasma at 
infinity, so we should expect our solution to radiate no flux. From ( 12.51) . it is clear that 



Regularity on the past horizon requires Q = 2K . 



the total flux radiated to infinity from the black hole is a constant proportional to K, 
so we conclude that K = and so Tt 2 ) makes no contribution to our stress tensor. 

• Unlike the past horizon, which is present only in maximal analytic extensions of black 
hole spacetimes, the future horizon is a genuine physical location present in any real- 
istic model of gravitational collapse. We thus expect physical quantities to be regular 
there. In particular, we require the stress tensor of the jammed CFT to be regular 
there; this implies that Q = 0, so that Tu\ does not contribute to the stress tensor 
either (incidentally, since we also have K = 0, this means that the stress tensor will 
be regular at the past horizon as well). 

We therefore conclude that only contribution to the stress tensor of the CFT state dual to the 
Schws droplet should come from f)2.9cl) . and so should only depend on the one function 0(r). 

In fact, we can go further and make claims about the behavior of 0(r). The free field 
results summarized earlier in this section should not all apply to our stress tensor, but we 
might draw certain universal behavior from them to make a guess at the behavior of O(r) 
(which we will then verify with our numerics). First, note that the negative energy density 
near the horizon present in free field theories is an indication of the highly non-classical 
nature of the field there; only far from the black hole does the stress tensor of the free 
fields become classical. But the jammed phase of the CFT is a highly non-classical state, 
as it consists of a plasma "halo" surrounding the black hole. In analogy with the free 
field theory results, we might therefore conjecture that the energy density of the strongly 
coupled CFT becomes negative near the horizon as well. In that case, we expand (T^) l t = 
— 3G(ro)/2 + Oir — r^) to see that a negative energy density near the horizon implies that 
is negative there. 

Far from the black hole, one can argue that due to the weak coupling between the 
black hole and the heat bath at infinity, the components of the stress tensor should fall off 
"rapidly". This statement can be quantified by following the logic of [1.1] : there, the authors 
find that the v.e.v. of a CFT dual to a d = 4 Schwarzschild droplet exhibits a 1/r 5 falloff. 



In order to explain this falloff, the authors invoke the results of 28], in which the linearized 
gravitational field created by a point-like source is calculated in the context of the Randall- 
Sundrum single braneworld model. By using the behavior of the gravitational field far from 



10 



the source, Einstein's equations 

R„u - 2 R9„u = MG N {T^) (2.10) 

can be used to obtain the expected large-r behavior of the stress tensor. For the case of 
a d = 4 dimensional spacetime, the predicted falloff of the components (T^ v ) goes like 1/r 5 , 
consistent with the results of |lj. One can generalize this argument to show that we ex- 
pect a l/r d+1 falloff in arbitrary dimension. This implies that for our solution, we should 
have 0(r) ~ 1/r 6 and G(r) ~ 1/r 3 at large r. 



B. A Little Spin 

Finally, let us consider the effect of giving the boundary black hole a nonzero angular 
momentum, thus "spinning" the droplet. Unfortunately, an analysis similar to that per- 
formed for the nonspinning case is not illuminating, because even with the conservation 
equations (j2.ip . there are too many independent functions in the stress tensor for a gen- 
eral solution to be tractable. Nevertheless, we can make some qualitative claims based on 
physical arguments. 

First, we clearly expect there to be angular flux terms as the CFT plasma is dragged 
around the spinning black hole. This flux density should be maximal near the horizon, and 
decrease monotonically away from the black hole. Similarly, the fact that the CFT plasma is 
being forced to rotate around the black hole leads us to expect that a centrifugal barrier forms 
around the black hole. Indeed, by considering timelike geodesies in the equatorial plane of 
the five- dimensional equal-angular-momentum Myers- Perry black hole, we can check that the 
centrifugal barrier of the effective radial potential grows and moves away from the horizon as 
the spin of the black hole is increased. This barrier may act to confine the plasma near the 
black hole horizon, effectively acting as a box around the black hole. We might expect that 
this effect would present itself as an increase in the radial and tangential pressures near the 
horizon. We should similarly expect an increase in the magnitude of the (negative) energy 
density near the black hole, as a buildup of negative-energy modes forms in the centrifugal 
box. This physical reasoning leads us to conjecture that all components of the stress tensor 
should increase in magnitude as the angular momentum of the black hole is increased. 
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C. Extremal Horizons 

Our family of solution black holes can be taken all the way to extremality. This limit is 
particularly interesting, as the causal structure of an extremal black hole is qualitatively very 
different from its non-extremal relatives. One might therefore imagine that the stress tensor 
of matter fields on extremal black hole spacetimes exhibits qualitatively different behavior 
as well. In particular, the fact that the horizon of an extremal black hole is a Cauchy horizon 



would naively lead one to believe that stress tensors should generically be singu 



ar on the 



future horizon of extremal black holes. This question has been addressed by |29l-l32j] . who 
studied the regularity of the v.e.v. of the stress tensor outside an extremal dimensionally 
reduced two-dimensional Reissner-Nordstrom (RN) black hole. In short, a static stress tensor 
(i.e. one sharing the same isometries as the background spacetime) exhibits a mild singularity 
on the future horizon, whereas the stress tensor of a massive scalar field propagating on an 
extremal two-dimensional RN geometry forming via gravitational collapse is regular thanks 



to the presence of decaying (but nonzero) flux terms. As shown in [33|], this subtle issue 
disappears in four dimensions, as the stress tensor of a scalar field in the zero-temperature 
vacuum state on the full four-dimensional RN spacetime is regular everywhere without the 
need for flux terms. 

To our knowledge, the stress tensor of a field theory on a rotating extremal black hole 
background has not yet been studied even in free field theory. The four-dimensional RN 
results might lead us to expect that the stress tensor on extremal rotating black hole space- 
times should similarly be regular on the future horizon, though a key difference between the 
RN and rotating case is the inherent presence of matter in the RN spacetime. In any case, 
the fundamental similarity between the causal structures of extremal charged and rotating 
black holes makes it quite plausible to expect that the CFT stress tensor dual to our rotating 
droplet will be regular on the future horizon even in the extremal limit. 
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III. CONSTRUCTING A SPINNING DROPLET 

A. The DeTurck Method 

The standard approach in AdS/CFT is to solve the vacuum Einstein field equations with 
negative cosmological constant 4 

Kb = Rab~ ~ Rgab + ^9ab = (3.1) 

subject to certain boundary conditions. The subleading behavior of the metric g a b near the 
conformal boundary then contains information about the stress tensor of the dual CFT. As 
is well known, Einstein's equations do not have a well defined character unless a gauge choice 
is made. Because the solutions we are searching for are stationary, one can hope to choose 
a clever gauge where the equations (13.11) . or a deformation thereof, are manifestly elliptic. 
This is exactly what the DeTurck trick does for us |l|, 



35] . In short, one modifies the 



equations (13. ip by introducing a new vector £ a (called the DeTurck vector): 

E» ee E ab - V (B &) = 0, C = 9 bc K ~ r y , (3-2) 

where f£ c is the Levi-Civita connection of some reference metric g. Equation (I3.2p is called 
the Einstein-DeTurck or harmonic Einstein equation. One can show that for stationary 
solutions with Killing horizons, the above choice of the DeTurck vector renders the Einstein- 
DeTurck equation elliptic. Solving the Einstein-deTurck equations then reduces to solving 
a boundary-value problem. Note, however, that solutions to the Einstein-DeTurck equation 
are not solutions to the ordinary Einstein equations unless £ a = 0. To get around this 
problem, one can show that the quantity $ = £ a £ a must take its maximum value on the 
boundaries of the domain of integration; thus if the reference metric g is chosen with the 
same boundary conditions as the metric g, $ is zero on all boundaries of the integration 
domain, and must therefore be zero everywhere within as well (l|. Then £ a is zero as well, 
and a solution to the Einstein-DeTurck equation is also a solution to the Einstein equations. 
Indeed, the claim that the Einstein-DeTurck equation automatically takes care of gauge- 
fixing is justified by thinking of the condition £ a = as just a choice of gauge, which is a 
generalized harmonic gauge of the form Ax a = T% c g bc . 



4 Here and below, we will use lower-case Latin letters to denote bulk indices, and lower-case Greek letters 
to denote boundary indices. 
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B. Droplet Ansatz 



As mentioned at the end of Section [J, the boundary metric we want to impose is the 
equal-angular-momentum d — 5 Myers- Perry solution 161 ]: 



ds 2 



MP 



gir) 

~h(r) 



dt 2 + 



dr 2 

g(r) 



+ r z 



h(r) (# + A(i) + n(r)dt) 2 + rfS 



h{r) = 

Q{r) = 



( r - r o)( r -P ( r + r o)) 
(1 - /3 2 )r 4 

r 4 - /3 2 (r 4 - rl) 

(1 - /3 2 )r 4 ' 

/3r 3 



-4\ ' 



(3.3a) 

(3.3b) 
(3.3c) 
(3.3d) 



< -/3 2 (r 4 -*o, 

where dT, 2 is the Fubini-Study metric on CP l , An.) is the Kahler potential of CP 1 , and /3 = 
tqQh, with f2# the angular velocity of the horizon. Explicit forms for An.) and dY\ are given 
by the expressions 



An 



1 



(i) 



cos 6 d<f), dill = -(d6 2 + sin 2 6 d<j) 2 ) 



(3.4) 



For (3 = 0, the metric ( I3.3P reduces to d = 5 Schwarzschild ( 12. 2D with the S" 3 written 
as a Hopf fibration of S" 1 over CP 1 ; setting /3 2 = /3 2 xt = 1/2 yields the extremal solution. 
The metric is cohomogeneity one owing to its U(2) symmetry, which implies that our bulk 
solution will be cohomogeneity two. 

Let us compactify the range of the radial coordinate r to yield a metric more amenable 
to our numerical approach. We also rescale t — > rot and exploit the coordinate freedom 
ip — > ip + Xt to shift Q so that fi(ro) = 0, fi(oo) = —Qh (this will simplify the boundary 
conditions at the horizon later on). Defining a new radial coordinate y = 1 — (r /r) 2 , we 
find 



ds 2 



MP 



-y 



hl{y) dt 2 + 



(3 2 )dy 2 



Mv) 



1-3/ 



4y/ii(y)(l-y) 3 

h 2 (y) 



1-/3 



:M ^ + A (1) -/3(l-/3 2 ) 2/(2 y) ^\ ^ 



h2(y) 



dt + dZi 



(3.5a) 



h 1 (y) = l-/3 2 (2-y), 
h 2 {y) = l-^y(2-y). 



(3.5b) 
(3.5c) 
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Our ansatz will need to have (J3.5P as its conformal boundary. 

Next, recall that we are considering black droplet solutions with T^ = 0; this implies 
that far into the bulk, the solution should approach the near-horizon geometry of pure AdS@ 
in Poincare coordinates. Poincare AdSg in standard coordinates can be written as 

ds 2 Poin = ^ [-dt 2 + dz 2 + dr 2 + r 2 dnl] , (3.6) 

where £ is the AdS 6 length scale, related to the cosmological constant by £ 2 = — 10/A 6 . 
Changing variables to (x, y) defined by 



the Poincare metric becomes 

f 2 
d* 2 - 

""'Poin 



1-X 



2^2 



1-X 2 



[l-y)dt 2 + 



x\ 



'2-x 2 



(3.7) 



dy 2 



+ 



4dx 2 



4(1 -y) 2 2-x 2 

The near- horizon region corresponds to taking y close to 1. 
We are now ready to write down an ansatz. Consider 

o2n „a r u t„.\ Ady 2 



+ x 2 {2-x 2 )dVtl 



(3.8) 



ds 2 



[i-v) f ..rrhi(y) Jl2 



'1-x 



2A2 



~ yT hZ) dt + lyh 1 (y)(l-y)* 



+ 



[l-y)(2-x 2 ) 



x(l - X 2 ) . 

dx + F -±- ^ dy 

2(1-1/) 



x 2 (2-x 2 ) 



h 2 (x,y)c(d^ + A {1) -G y ^ ' y) 



dt) +SdY>l 



where 



h 2 (x,y) = 1 + x 2 



h 2 (y) 



(3.9) 



(3.10) 



1-/3 2 

T, A, B, C, S, F, and G are all functions of x and y, and the coordinate range is the 
rectangle (x,y) £ (0,1) x (0,1). For the non-extremal case j3 ^ l/y/2, our boundary 
conditions are as follows: 



At the conformal boundary x — Y 1, we take T = A/(l — (3 2 ) 
1, G = (3(1 — (3 2 ), and F — 0. Then the metric approaches 

' J ' 1-3/ 



C = 5 = r 2 ,/^, 5 



rfs 2 ^ 



X 



<ix z 



,7,2 

4 ^ 2 GS MP 



(3.11) 



as desired. 
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Near the extremal Poincare horizon y — > 1, we take T = A /(l — (3 2 ) = B = C = S = 
1, G = f3(l — (3 2 ), and F — 0. Changing to a new coordinate ^ — > ip + /3t removes 
the cfa cross-terms, and we recover the y — > 1 limit of (13. 8p . 

At the horizon y — > 0, we require regularity of the metric, which imposes a rela- 
tionship between T and A. In particular, we must have T/A = (ro/t/(l — 2/3 2 )) 2 , 
where k, — (1 — 2/3 2 )/(roA/l — /3 2 ) is the surface gravity of the boundary black hole. 
Thus we find T/A = 1/(1 — /3 2 ). In addition, expanding the equations of motion and 
the condition £ a = order-by-order near y = gives Robin conditions on d y Mi\ y= o 
for M t = {T, A, B,C,S,F,G}. 

At the fixed point of the U(2) isometry x = 0, we require regularity as well. Again, this 
implies that all metric functions must be smooth functions of x 2 , giving the Neumann 
conditions d x Mi\ x= o = 0. In addition, note that while the forms x 2 (dip + A^) and xdx 
are regular at x — 0, in general the metric components 



2 B dx 2 + 2 x 1 



h 2 (x,y)C(dilj + A {1) ) 2 + SdJ:l\ (3.12) 



are not. In order to make the above expression regular, we must impose that B 
C = S at x = 0, so that 



2 B dx 2 + 2 x 2 



h 2 (x,y)C(dtlj + A {l) ) 2 + Sd£ 2 ->■ 2S(tfe 2 + a; 2 ^ 2 .), (3.13) 



which is manifestly regular there. 

Having understood the boundary conditions on each side of the computational rectangle, 
we choose the reference metric. Noting from H 3 . 5 1) that the boundary black hole size vq 
enters only as a conformal factor, we can choose without loss of generality to set ro = £. 
Then the boundary conditions become consistent, and we take the reference metric to be 
the ansatz flU} with T = A/(l - (3 2 ) = B = C = S = 1, G = (3(1 - /3 2 ), and F = 0. 

The extremal case ((3 = (3 cxt = 1/ ' \/2) takes a little more care. Though our ansatz (13. 9p 
still applies, one finds that as extremality is approached the metric function F attains 
a very large value along y = near x = 0, leading to large numerical errors. To 
alleviate the problem, we rewrite (I3.9P in terms of a new metric function F = yF. 
Then the boundary conditions at the conformal boundary, fixed point of the U(2) isom- 
etry, and extremal Poincare horizon remain the same as those discussed above, while at 
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the extremal droplet horizon we require F = 0. In addition, converting to Eddington- 
Finkelstein coordinates and requiring regularity at the extremal droplet horizon imposes 
the usual regularity condition A(x,0) = (1 — j3^ xt )T(x,0) = T(x,0)/2, as well as the con- 
dition G(x, 0) = const. = /3 ex t(l — /3g X t) = 1/V&- Boundary conditions on the other metric 
functions are obtained by expanding the equations of motion order-by-order near y = and 
requiring that they vanish to leading order. 

C. Numerics 

We approximate the Einstein-DeTurck equation ( 13.21) with the ansatz (J3.9P using pseu- 
dospectral collocation methods on a Chebyshev grid. The resulting nonlinear algebraic 
equations are solved using a damped Newton-Raphson method. We monitor the damping 
with |-E^|; that is, we ensure that each iteration of the Newton's method decrease the mag- 
nitude of E^ b . Due to the maximum principle that $ = £ a £ a obeys, we can monitor the 
error in our solutions by monitoring the maximum value of $. In Figure [21 we plot the 
maximum value of $ as a function of the number N of grid points for three different values 
of the rotation parameter /5, showing the expected exponential decrease. Note that we will 
find later that the accuracy of the extracted of the stress tensor depends strongly on the 
number of grid points used along the re-direction, so for all results shown in this paper, we 
use an N x x N y = 81 x 41 grid. 

D. Extraction of the Stress Tensor 

In order to extract the boundary stress tensor from our solutions, we follow the approach 



used in [l3|. In short, we first expand the metric functions Mj as a power series in x off of 



the boundary at x = 1: 

M,(x, y) = J2 Mt\v)(l - *T + Ml - *) E Mt\v)(l - xY, (3-14) 

ra=0 n=0 

where the logarithmic terms only appear for (3 ^ and at no lower order than n = 5. 
Inserting these expansions into (I3.2p and in addition imposing £ a = allows us to solve for 
the coefficients of the power series order- by-order in (1 — x). The expansion is unique up 
to n = 4 for the functions {T, A, C, S, F, G} and up to n = 5 for the function B. We then 
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FIG. 2: The maximum value of the square of the DeTurck vector $ as a function of N for (3 = 
(circles), 0.6 (diamonds), and (3 ext ~ 0.707 (asterisks). 



change to Feffer man- Graham coordinates (z, y) using an expansion of the form 



i-x 2 = v/r^+^rx(")(y)^, 

n=1 



(3.15a) 
(3.15b) 



?1=1 



where the O(z) term in the expansion for x fixes the conformal frame, and we have neglected 
potential logarithms as they do not appear up to the order we need. Requiring that g$% = 
£ 2 /z 2 and g^ y = order-by-order in z yields a set of algebraic equations that can be solved 
for the coefficients X^ n '(y), Y( n \y). Then we find that the metric takes the form 

p2 

|2 



ds 2 = ^ [dz 2 + £~ 2 ds 2 MP + th^ dx»dx v + 0(zf] 



(3.16) 



where h^ v is only a function of the boundary coordinates x M . According to the prescription 
of 36j, the renormalized v.e.v. of the boundary stress tensor is then 

5£ 4 



{T,u) 



■I « ^(6) h W 
107TGJY 



(3.17) 



The general expressions for the coefficients M\ , M\ n , X^ n \ Y^ 1 and for the stress tensor 
are too cumbersome to reproduce in this paper, but they become tractable in the nonrotating 
case (3 = 0. In Appendix [Aj we give the (3 = expressions for M\ , X^ n \ and Y^ n ' up to 
the order we need in order to extract the stress tensor and check its tracelessness and 
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transversality. As expected from our arguments in Section III Al the stress tensor is given 
precisely by the expression (I2.9cl) . with 5 

5 £4 (1 _ y) 5/2 [M(1 _ 2 y)a 5 (y) + (1 - j/)(176 - 1712y + 2368y 2 + 14^0/))] 



©(</) 



IGttG^ 672(1 - 2y) 



(3.18a) 
^ , 5l 4 y(l - y) (88 - 38% + 296y 2 + 7a 5 (y)) 
G(2/) = To^f 2^4 ' (3 - 18b) 

where in an obvious abuse of notation we have written the functions and G in terms of 
the compactified radial coordinate y rather than the original coordinate r. We will postpone 
a discussion of this result to Section HVl 

Now let us return to the general /3 7^ case and discuss how to extract the stress tensor 
from our numerical data. The boundary stress tensor will take the form 

(T^)dx»d u = T tt dt 2 + T rr dr 2 + T^ (dip + A {1) ) 2 + 2T H (dip + A {1) ) dt + T ES rfS 2 . (3.19) 

There are thus five functions to solve for, which can be expressed in terms of five of the 
functions M- which we will label as m\ = {t 5 , a 5 , c 5 , g 5 , s 5 }. In principle, one could extract 
these functions by taking five ^-derivatives of (13.141) at x = 1. However, in practice taking 
such high derivatives numerically gives very poorly behaved results. Instead, we follow a 

(n) 

slightly different approach: since we know the coefficients M\ for n < 4 analytically, we 
can write 

(1 — x) 4 

With the numerically computed values for Mi(x,y), the left-hand side of (I3.20p is known 

(5) 
numerically, so we perform a fit for the right-hand side to obtain the coefficients m\ . 

In order to quantify the accuracy of our extracted stress tensor, we proceed as follows. 

Although the equations of motion don't give unique expressions for the coefficients m\ , 

they do provide two algebraic relations between these coefficients and a' 5 6 . Using one of 

these relations, we can express £5 in terms of 05, C5, and g$ in an algebraic expression of the 

form 

h(y) = Tr ct (V, 05(2/), c 5 (y), s 5 (y)) . (3.21) 



5 Note that we drop the tilde on y, since y and y agree on the boundary. 

6 These relations enforce tracelessness and transversality of the stress tensor, i.e. conservation. 
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FIG. 3: The relative error A introduced in extracting the coefficients m] for different values 
of j3. Circles, squares, diamonds, triangles, inverted triangles, and asterisks correspond to /3 = 
0, 0.5, 0.6, 0.65, 0.69, and (3 ext respectively. Note that the error remains < 0.4% for almost all f3 
and y. 

Using the fitted values for the mf , we can therefore use the quantity 



A(y) 



h(y) 



(3.22) 



as a measure of the error introduced in extracting the m\ (physically, one can think of A 
as a measure of how much the numerically extracted stress tensor fails to be traceless) . In 
Figure [31 we plot A for the full range of /3. We note that the relative error in our extraction 
is < 0.4% for almost all /3 and y; only for y near 1 does A become appreciable. We emphasize 

(6) 

that this error is solely due to the difficulty in extracting the m\ from the numerical data, 
and is not a measure of the accuracy of our numerical solutions themselves. 

In Figure HI we plot the components of the stress tensor for various values of (5. Recall 
that to obtain the metric (I3.5p . we shifted the angular coordinate ip to make g^t = at 
the horizon. As a result, the near-horizon behavior shown in Figure H] physically represents 
the stress tensor as measured by an observer co-rotating with the horizon. In addition, in 
Figure El we plot the scalar invariant (T^ u ) (T^ u ). 
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FIG. 4: The components h^ v = 16irG N (T^ v ) /5£ 2 for the same values of /3 as Figure Note 
that h 1 ^ = for /3 = 0. Note that the range of the coordinate y shown corresponds to the 
range tq < r < 1.2 ro in the original radial coordinate. 
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1. Regularity at Horizon 



From the discussion in Section Til Al and the regularity of the functions (j3.18p . it is clear 
that in the nonrotating case, the stress tensor of the CFT is regular on the future horizon 
(and in fact, it is regular on the past horizon as well since K = Q = 0). One can similarly 
examine the behavior of the stress tensor in the rotating case by passing to local ingoing 
Eddington-Finkelstein coordinates. In the non-extremal case, we change coordinates to 



1 - /3 2 dy 



dt = dv- V "I ^; (3.23) 

1-2/3 2 2y V ; 

then from f 1 3 . X 9 1) . one finds that the stress tensor is regular on the future horizon only if all the 
components (T^ u ) are finite and satisfy (T* t ) = (T y y ) on the horizon. We can immediately 
see from Figure H] that these conditions are satisfied, so the stress tensor is regular on the 
non-extremal horizon. 

Again, the extremal case requires more care. In this case, local ingoing Eddington- 
Finkelstein coordinates are given by 



dt = dv = ( — - H — I dy, 

dip = dip 1 . 

2y 



(3.24a) 
(3.24b) 



Then one again finds that the stress tensor is regular on the future horizon if it obeys the 
same conditions as in the non-extremal case, and in addition obeys 



V2(T t lp ) + 2d y {(T t t )-(Ty y ))=0 
22 



(3.25) 



on the horizon. It is easy to check that this condition is satisfied to within the expected ~ 1 % 
accuracy. 

IV. RESULTS AND DISCUSSION 

Our numerical results exhibit the behavior expected based on our quantitative and heuris- 
tic reasoning in Section [III Let us begin with the nonrotating case /3 = 0. The expansion of 
the equations of motion off of the boundary gave us a stress tensor of the form f )2.9cp with 
the functions and G given in f)3.18p . which is precisely the form expected from the discus- 
sion in Section III At In particular, we find that the stress tensor is traceless, contains no flux 
terms, and is regular at the future horizon. Once the numerically extracted functions mf(y) 
are inserted, we find that the energy density is negative everywhere. The magnitude of 
the components of the stress tensor increases as the rotation parameter is increased, but 
otherwise the qualitative behavior of the stress tensor doesn't change when /3 ^ except for 
the introduction of an angular flux term T t ^. 

As shown in Figure [3j the error in our extracted stress tensor becomes relatively large 
near y — 1, making it difficult to reliably extract the falloff of the stress tensor components 
at large r. However, studying the behavior of the stress tensor components for which A 



is smallest indicates a fa 
braneworld arguments of 



loff consistent with the (1 — y) 3 = (ro/r) 6 expected from the 

The following physical pictures emerges. As expected, a strongly coupled large- N CFT 
in a jammed phase forms a halo of negative-energy plasma around a black hole. The stress 
tensor of the plasma falls off rapidly far from the black hole, indicating that the plasma is 
well-localized around the black hole. The total energy of the plasma is negative, indicating 
its highly non-classical nature. When the black hole is made to rotate, the plasma is dragged 
along with the black hole. The plasma is trapped by the centrifugal barrier created by this 
rotation, causing the energy density and pressures to increase in magnitude. 

The jammed phase is particularly interesting from the point of view of the CFT because 
it has no analog in a free field theory: a similar static plasma localized around the black 
hole could not exist in a noninteracting theory, as it would quickly fall into the black hole. 
The jammed phase is thus an effect of the strongly coupled nature of the CFT. 

One might therefore wonder how our droplet phase compares to the corresponding black 
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funnel phases (which do admit analogs in free field theory). Black funnels with T^ = 
must have a horizon that asymptotes to the extremal Poincare AdS deep into the bulk; if we 
require that the bulk horizon also join smoothly to the (non-extremal) boundary horizon, 
the bulk horizon must have both a non-constant surface gravity and a non-constant angular 
velocity, and will therefore not be a Killing horizon. We say that such funnel solutions will 
"flow" in the sense that there exists some notion of a horizon velocity which moves either 
from the boundary black hole to infinity or vice versa. For non-extremal horizons, these 
flowing funnels were constructed in 14|, but construction of solutions with asymptotically 
extremal horizons (as is necessary to have the boundary CFT be at zero temperature at 
infinity) is more difficult. 

In fact, it would seem that such funnel solutions might not even be relevant: according 
to the discussion in Section [J the conjectured jamming transition in the CFT should be 
parametrized by the parameter RT^, with small RT^ favoring the jammed phase. If T^ = 0, 
we might expect based on this argument that the jammed phase of the CFT will dominate 
the thermodynamic ensemble for any size of the boundary black hole. 

However, recall that one basis for this conjectured phase transition was the potential for 



Gregory-Laflamme-type instabilies of the bulk horizon 37|: if one starts with a black funnel 
and lowers the temperature T^, the asymptotically planar bulk horizon will sink deeper into 
the bulk (cf. Figured]). The naive expectation is that this will narrow the neck of the black 
funnel and leave it unstable to Gregory-Laflamme instabilities, causing it to collapse into 
a black droplet. Taking T^ — > 0, this argument would lead one to expect that the droplet 
solution would be the only thermodynamically preferred one 



This potential instability was studied to some extent in [13], which considered funnels 
dual to a CFT on an asymptoticaly flat black hole spacetime. The size of the boundary 
black hole was decreased (while keeping the temperature T^ fixed) to explore whether such 
Greogory-Laflamme instabilities would occur. Surprisingly, as the size of the boundary black 
hole decreased, the neck of the funnel remained wide enough to prevent the formation of 
such instabilities. Though far from exhaustive, this result implies that it is possible for there 
to be stable funnels with T^ = which could compete with the droplet solutions. A sketch 
of what such a solution might look like is shown in Figure [6] (in fact, in [l2[ such funnels 
were explicitly constructed in d = 2, though the question of their stability is less interesting 
in that case, since the droplet solutions cannot exist in three bulk dimensions). Clearly, the 
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FIG. 6: A sketch of a black funnel solution with T^ = 0. 

construction of such solutions would be of much interest. We leave this exploration to future 
work. 

Some words on the extremal limit are also in order. Though this limit is discontinuous 
from the point of view of the background geometry, the stress tensor changes smoothly 
as (5 approaches /3 ext . In particular, the stress tensor remains regular on the future horizon, 



consistent with the result of 33( for the four-dimensional RN spacetime. A natural question 
is then whether there is some other mechanism that prevents the spacetime from continuing 
past the Cauchy horizon, analogous to the mass-inflation singularity of black holes with 



inner horizons J38|; perhaps the inclusion of perturbative corrections at higher order in 1/N 
would resolve the issue. 

Indeed, our results may also have implications for the mass-inflation singularity. Ref. \±[ 
computed the stress tensor on a rotating BTZ black hole background and found that it 
diverged at the inner horizon. This result should be expected in any black hole spacetime 
with an inner (Cauchy) horizon to preserve the causal structure of the spacetime. Since 
the spacetime between the inner and outer horizons is dynamical, the numerical approach 
used in this paper would not be applicable to studying the event horizon in that region. 
However, obtaining the stress tensor of a field theory on a black hole background all the 
way to the inner horizon would certainly give some insight into the behavior of the mass 
inflation singularity there. 
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Appendix A: Near-Boundary Expansion for the Nonrotating Droplet 

In this appendix, we list the expansions of the metric functions near the conformal bound- 
ary x = 1 for the nonrotating ((3 = 0) droplet, as well as the form of the Fefferman- Graham 
expansion (I3.16p . For the metric functions, we have 



T{x,y) = l + {l-y) (-6(l-x) 2 + 6(l-x) 3 + - y ' 1 - ! 



:i-x) 



14 

+ t 5 (y)(l-xf + 0(l-xf, (Ala) 



A(x, y) = i + (i-y)\ -6(1 - xY + 6(1 - xf + 



3 , 67-296?/%, ^, 



14 



1-x) 



+ a 5 (y)(l-xf + 0(l-xf, (Alb) 



S(x, y) = C(x, y) = 1 + (1 - y) ( 2(1 - xf - 2(1 - xf + 5 ( 19 + 8 ^) (1 - x)^j 

+ s 5 (y)(l-xf + 0(l-xf, (Ale) 



B(x, y) = l + -(l-y)(5 + y) ((1 - xf - 2(1 - xf) + 0(1 
F(x,y) = -(l-y)(l-xf 



— x) 



s 211-192?/, ,n /n/ ,o 

1 + (1 - x) + ^^ (1 - xf + 0(1 - xf 



G{x,y) = 0. 



(Aid) 
(Ale) 
(Alf) 



The coefficients t$(y), a 5 (y), and Ss(y) cannot be determined uniquely, but must satisfy the 
relationships 

7(1 - 2y)a B (y) - (1 - y)(Q5Qy - 9Wy 2 - Uya' 5 (y)) 



h{y) 
Sb{y) 



7(1 - 2y) 
14(1 - 2y)a 5 {y) + (1 - y)(8(55 - 265y + 266y 2 ) + 14y^(j/)) 



21(1 -2y) 

For the transformation to Fefferman- Graham coordinates, we have 



1-x = z^\ - y 



1- 



\_-yf_ 3(l-^(3 + 6y + 19y 2 ) 

4 + 224 + l j 



(A2a) 
(A2b) 

(A3a) 
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y = y 



-2,2 , (i-yf(i-Sy) z4 



l + (l-y) 2 z 2 + ± ^ ^5 4 

(A3b) 



1 - //) '( 1 -■- 19y - 44y 2 ) ^ | 7 
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